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Abstract 

The available upward-going muon data from the Kamiokande, Baksan, 
MACRO, 1MB, and SuperKamiokande experiments are reviewed and com- 
bined. Bounds on the neutrino mass and mixing parameters are derived for 
oscillations in two and three flavors. These bounds are not in significant con- 
flict with the oscillation solution to the atmospheric neutrino flavor anomaly 
observed in the sub-GeV and multi-GeV energy range. The combination of 
all the available atmospheric data tends to favor the ^ Ve channel with 
respect to the <^ Uj- channel, and to disfavor the threefold maximal mixing 
scenario. 

PACS number(s): 96.40.Tv, 14.60.Ef, 25.30.Pt, 14.60.Pq 
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I. INTRODUCTION 



The observation of upward-going muons in underground detectors has long since been 
considered |l| as an effective tool to study the neutrino component expected in atmospheric 
showers 0-^. Early anomalies in the muon event samples detected by the first two dedicated 
experiments soon sparked investigations of possible neutrino oscillation effects 0J^. 

The upward-going muon total rates measured in recent underground experiments appear 
to be approximately in agreement with the expectations, thus disfavoring large neutrino 
oscillation effects. This contrasts with the observed anomaly in the electron and muon 
contents of (semi) contained events induced by atmospheric neutrinos, that might signal 
flavor oscillation effects (see ^ for a recent review). The conflict between these two data 
sets is usually resolved by appealing to the large, noncancelable uncertainties affecting the 
calculation of the absolute muon rates [[10|-[T^ .[| In terms of neutrino oscillations, this implies 



that there are values of the oscillation parameters that explain the "anomalous" muon-to- 
electron flavor ratio without really conflicting with the "regular" upward muon data |[r3|-[T5|. 
In this work we confirm and elucidate such point of view by performing an updated analysis 



of the upward-going muon data from five experiments: Kamiokande p6|-[T9||, Baksan p0|-p2|, 
MACRO 1MB 0,0], and SuperKamiokande [0,^ . 

The paper is structured as follows. In Sec. II the available data on upward-going muons 
are reviewed, compared with each other and with the theoretical expectations (in the absence 
of oscillations), and then combined. In Sec. Ill the combined data on the muon flux angular 
distribution are used to derive bounds on the two-flavor and three-flavor neutrino oscillation 
parameters. These bounds are compatible with the neutrino oscillation solution to the 
atmospheric neutrino anomaly. It is shown that the inclusion of upward-going muon data 
in fits tends to favor the z/^ ^ oscillation channel and to disfavor the threefold maximal 



^ We add that the experimental systematics might also be more significant than usually assumed, 
see Sec. II B. 
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mixing scenario. In Sec. IV we draw the conclusions of our work. This study is part of a 
vast research program aiming to analyze the world neutrino oscillation data, including the 



results of solar ||2^, atmospheric, [p9|-|3T[] and laboratory (accelerator and reactor) p9|j32 
neutrino experiments. 



II. EXPERIMENTAL DATA AND STANDARD EXPECTATIONS 

The experimental data on upward-going muons are usually given as distributions of 
either muon events or muon fluxes as a function of the zenith angle 9, with cos0 = 
(cos0 = —1) corresponding to horizontal (vertical) muons. The distribution of muon events 
depends upon technical specifications of the detector such as the efficiency and the geometric 
acceptance, that are often unpublished. The distribution of muon fluxes (when available) 
is instead deconvoluted for these effects, and thus allows a more direct comparison with the 
theoretical expectations. We describe in Sec. II A the more recent available data (either 
muon events or muon fluxes) and in Sec. II B the specific combination of data used in our 
oscillation analysis. 



A. Data from different experiments 

In the last few years, data on upward going muons have been collected by five ex- 



periments: Kamiokande lT|-|Ti, Baksan MACRO 1MB |2|,|2|, and Su- 

per Kamiokande p7|j28[| . The 1MB experiment (completed) was performed in three distinct 
stages (phases 1, 2, and 3); we will comment both on the global results (IMB-1,2,3) p6[ and 
on the results from the first two phases (IMB-1,2) [|25|. Data from the pioneering NUSEX 
experiment (eight events only) will not be considered. Since some experimental results 
have been reported only in conference proceedings, we compile the most recent data on 
upward-going muons for the benefit of the reader. The compilation is reported in Table | 
and in Fig. |1|, that we comment in parallel. 
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Table | reports, for each experiment, the relevant detector characteristics and the ex- 
perimental muon data, compared with the corresponding expectations in ten angular bins. 
The theoretical muon fluxes refer to our calculation^ with the following inputs: Bartol neu- 
trino fluxes p6| (neglecting small, site-dependent geomagnetic effects), version "G" of the 
Martin-Stirling- Roberts [MRS(G)] structure functions [^] (i.e., the default choice of the 
CERN Parton Density Function Library [^,^), and Lohmann-Kopp-Voss (LKV) muon 
energy losses in standard rock [|^. The use of the same input for the calculation of the 
expected muon fluxes will enable us to perform a uniform comparison with all the experi- 
ments.0 Notice that, for flxed inputs, the theoretical muon spectrum depends only on the 
energy threshold; therefore, the Baksan and MACRO experiments (both with > 1 GeV) 
share the same expectations. All other experiments have higher thresholds^ and correspond- 
ingly lower muon fluxes. Concerning the experimental data, we quote only the statistical 
errors, since the systematic errors in each bin are usually unpublished. The only notable 
exception is the MACRO analysis p3| , p4| , where the systematic uncertainties are estimated 



to range from a few percent near the vertical to about 20% near the horizontal direction. 
We also report (last row in Table ^ the total muon flux, deflned as the sum of the individual 
fluxes in each bin times the bin width, with (statistical) errors added in quadrature. 

Figure |l| displays, in graphical form, the same information as Table H The data are shown 
as dots with (statistical) error bars. Our calculations (in the absence of oscillations) are 
represented by the solid histograms. In addition, we report as dashed histograms the original 



^ A general description of the theoretical estimate of muon fluxes can be found, e.g., in [34, 3^]. 

^It should be stressed, however, that alternative input choices for the neutrino fluxes or for the 
structure functions may induce variations as large as ~ 30% in the normalization of the expected 
muon fluxes. This uncertainty will be included in the fits. 

^In particular, the 1MB threshold corresponds to ~ 1.8 GeV for phases 1 and 2, to 1 GeV for 
phase 3, and to an effective value ~ 1.4 GeV for the total data sample (IMB-1,2,3). 
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MonteCarlo simulations of the various collaborations. The ingredients of these simulations 



were: (1) Kamiokande [0,|T3]: Volkova neutrino fluxes [O], Eichten-Hinchliffe-Lane-Quigg 



(EHLQ) structure functions p2|,|38|, and Lohmann-Kopp-Voss (LKV) muon energy losses 



P0| ; (2) Baksan [pT] , p2| : Bartol neutrino fluxes [^], Morfin-Tung (MT) structure functions 



and LKV muon energy losses; (3) MACRO [^|: Bartol neutrino fluxes, MT structure 



functions, and LKV muon energy losses; (4) 1MB ^,|6|: Volkova neutrino fluxes, EHLQ 
structure functions, and Bezrukov-Bugaev muon energy losses |0; (5) SuperKamiokande 
28| : Bartol neutrino fluxes, EHLQ structure functions, and LKV muon energy losses. The 



differences between our calculations and the MonteCarlo simulations are relatively small. 

The last panel in Fig. |l| shows the total number of muons that stopped in, or passed 
through, the 1MB detector (phases 1, 2, and 3 combined) .0 Unfortunately, it is difficult to 
calculate reliably the corresponding theoretical rates without detailed information about the 
1MB detector. The very definition of "passing" or "stopping" muon depends on unrepro- 
ducible geometric cuts applied on an event-by-event basis. Their approximation with (more 
manageable) energy cuts (see, e.g., [|n|) may bias the classification of events. Moreover, 
the absolute estimates of muon yields requires also the accurate knowledge of the detec- 
tor geometrical acceptance and of its energy-dependent efficiency, which are not published. 
Therefore, we are not surprised that, just for IMB-1,2,3 (last panel of Fig. |l]), our tentative 
calculation of the "passing" muon yield is not in good agreement with the published 1MB 



MonteCarlo simulation ||2^. For these reasons, we have not used the IMB-1,2,3 stopping 
and passing muon data in our analysis. We have instead used the muon flux distribution 
(published for IMB-1,2 p5[), that can be calculated more reliably and can also be compared 
with similar distributions from other experiments. 



The muon flux distribution for IMB-1,2,3 is, unfortunately, unpublished. 
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B. Combination of data 



As Fig. |l| shows, the angular distributions measured by Kamiokande, Baksan, MACRO, 
IMB-1,2, and SuperKamiokande are roughly in agreement with the theoretical expectations. 
However, the data also show some "bumps" or "dips" that seem too large to be statistical 
fluctuations but also too random to suggest new physics — they are likely to be the effect of 
unknown experimental systematics. Let us consider, e.g., the MACRO and Baksan experi- 
mental spectra, that, as observed in the previous subsection, are characterized by the same 
threshold {E^ > 1 GeV) and thus can be directly compared by superposition. Both the 
Kolmogorov-Smirnov and the tests give a <^ 1% probability that these two distributions 
are the same within the statistical errors. This suggests that the experimental systematics 
might be nonnegligible (in some bins, at least), as indeed confirmed by the MACRO error 
analysis [^. Unfortunately, no analogous error analysis has been reported by the Baksan 
Collaboration. 

The presence of relatively large "dips and bumps" in the individual experimental spectra 
of Fig. 1^ poses some practical problems in fitting the data with theoretical curves (that, 
conversely, are rather smooth even in the presence of neutrino oscillations). In fact, the 
oscillation fits to the individual spectra become dominated by the bins that exhibit the 
largest fluctuations, thus producing rather unstable and unreliable results. Therefore, it 
seems more reasonable to combine first the data from different experiments, in order to 
"average out" their fluctuations, and then to fit the neutrino oscillation parameters. 

We have combined the five muon flux spectra of Fig. |l] into a single spectrum as follows. 
First, the experimental distributions have been "rescaled" to a common 1 GeV threshold, by 
adding in each bin the theoretical contribution due to the change in threshold (no correction 
is necessary for MACRO and Baksan). The results are shown in Fig. 0(a), together with 
the theoretical expectation (solid histogram), now common to all experiments.^ Then 

^ Notice that the spread of the experimental data points in Fig. ^(a) is significant (a factor of about 
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the data have been combined in each bin, using only the known, statistical errors. The 
resulting, combined distribution is shown in Fig. ^(b) (dots with small error bars): it appears 
smoother than each of the individual parent distributions, and thus more appropriate for 
fits. A delicate point is the estimate of the (nonnegligible but unknown) systematic errors. 
The Particle Data Group recipe for scaling data errors as \Jx^ / {N-qy — 1) [^] would give 
scaling factors between 1 and 2 in the various bins. Similarly, the addition of plausible, 
(9(10%) systematic errors would give total errors ~ 1.5 times larger than the statistical 
errors. Guided by these estimates, we have chosen to rescale the statistical errors by a 
uniform factor 1.5 [large error bars in Fig. 0(b)] in order to obtain the "total" uncertainties. 
These "inflated" experimental errors have been used in all the following fits. 

We stress that any conceivable combination of the present muon data necessarily involves 
some approximations or subjective choices, both because spectra with different thresholds 
cannot be compared directly and because no experiment (excepted MACRO) reports the 
systematic uncertainties in each bin. For instance, the (somewhat arbitrary) choice of a 
common threshold for rescaling the experimental spectra has some influence on the sensi- 
tivity of the fits to low values of the neutrino mass differences. Therefore, great care should 
be taken in interpreting the results of any fit to the available muon data, including those 
presented in the following Section. Hopefully, in the next few years the SuperKamiokande 
experiment will collect much higher muon statistics and will carefully estimate the system- 
atic uncertainties of the flux distribution, so as to render such drawbacks gradually less 
significant in future analyses. 



two in the muon fluxes), suggesting again and independently that large and not well- understood 
systematics might affect (at least some of) the experimental results. 



7 



III. UPWARD-GOING MUONS AND NEUTRINO OSCILLATIONS 



In this section we report the results of our two-flavor and three-flavor oscillation analy- 
sis of the upward muon spectrum shown in Fig. 0(b). Since this combined spectrum is in 
reasonable agreement with the expectations in the absence of oscillations, exclusion zones 
can be drawn in the neutrino oscillation parameter space. This contrasts with the fla- 
vor anomaly observed in contained and semicontained atmospheric neutrino events, which 
suggests nonzero neutrino masses and mixings However, as we will see, positive and 
negative indications are compatible, within the experimental and theoretical uncertainties, 
in some zones of the parameter space where the neutrino mixing is "intermediate" (i.e., 
neither zero nor maximal). 

We use the three- flavor oscillation formalism defined in our previous works p9| , |3Tl , |32| . 
The first two of the three neutrino mass eigenstates (z/i, z/2, 1^3) are assumed to be essentially 
degenerate (mi ~ 7712), and separated from the third state z/3 by a square mass gap = 
Img — ml2\- The flavor content of z/3 is parametrized as 

z/3 = sin z/e + cos (sin ip + cos z/^) , 

where the mixing angles and ^p range between and tt/2. The cases (j) = 0, ip = 7r/2, 
and ip = correspond to pure two- flavor oscillations in the channels z/^ z/^, z/g z/^, and 
z/g z/,-, respectively (the latter case z/g ^ z/^ is uninteresting for atmospheric neutrinos and 
will not be specifically considered). As shown in [plj^ , it is useful to chart the parameter 



space through the logarithmic coordinates (m^, tan^ 0, tan^ ip). Our results will be presented 
in representative sections of this space. 

For each value of the neutrino mass and mixing parameters, we calculate the theoretical 
muon flux spectrum with the inputs described in the caption to Table |, including the 
matter oscillation effect for neutrino paths in the Earth p9| , |31[| . Then we evaluate the 
goodness-of-fit to the experimental distribution in Fig. ^(b) through a test. The 
quadratic form includes both the uncorrelated total experimental errors and the correlated 
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theoretical uncertainties of the absolute neutrino fluxes, cross sections, and muon energy 
losses in the rock. The Icr total theoretical error is conservatively assumed to be 30%, 
with a bin-by-bin correlation coefficient equal to +1. For both two-fiavor and three-flavor 
oscillations the minimum value of is searched, and CL. contours are drawn at the values 
of Ax^ = ~ Xmin appropriate to A^df = 2 and A^df = 3, respectively (see also [PT[]). 
We will combine in different ways the following data sets, in order to show their relative 
"weight" in driving the fit: upward-going muons (U) from Kamiokande, IMB-1,2, Baksan, 
MACRO, and SuperKamiokande; sub-GeV (S) yU-like and e-like events from Kamiokande, 
1MB, NUSEX, and Frejus [0; multi-GeV (M) yU-like and e-like events from Kamiokande 
46| , both binned (Mb) and unbinned (Mu) ||31|. The data sets S, Mu, and Mb are analyzed 



as in iSTII 



A. Two-flavor oscillations 

Figure ^ shows the results of our analysis of ^ Vr oscillations. In the first panel 
(U) only the upward-going muon data are fitted. The symmetry of the CL. contours with 
respect to the vertical axis = 7r/4 reflects the absence of matter effects in the z/^ ^ Vr 
channel ||3l[|. As expected, large neutrino mixing is excluded (for m? ^ 10~^). This negative 



indication, however, does not rule out the evidence for oscillations coming from sub-GeV and 
multi-GeV binned or unbinned data (S-|-Mb and S-I-Mu panels). In fact, the combination of 
all atmospheric data still favors the oscillation hypothesis, as evidenced in the plots labeled 
"S+Mb+U" (sub-GeV + multi-GeV binned + up-going n) and "S+Mu+U" (sub-GeV + 
multi-GeV unbinned + up-going /i). 

Notice that the addition of the upward-going muon data reinforces the exclusion of large 
^ mixing (required by the negative results of the NUSEX and Frejus experiments [pl[|), 
but does not alter dramatically the oscillation fit to sub-GeV and multi-GeV atmospheric 
neutrino data. In particular, for large vr? the "S+Mu+U" bounds on tan^ are only slightly 
narrower than those without upward /x data ("S+Mu" panel). This observation is relevant 
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for those phenomenological models in which the controversial zenithal anomaly of multi- 
GeV data is discarded, such as the "minimum sacrifice" scenario studied in P7| , ^ . We 
have verified that the solution to the neutrino oscillation evidence proposed in such scenario 
survives after the inclusion of the upward-going muon data, although with a slightly worse 

Figure § shows the results of our analysis of Ve ^ i^fj. oscillations. The panels are labeled 
as in Fig. Notice that the bounds from upward-going muon data in the ^ z/^ channel 
are considerably weaker than in the Vr channel, and basically exclude a relatively 

small zone at 90% CL. for large mixing and im? ^ 0.06-0.3 eV^ (no constraints appear at 
99% C.L.). In fact, in the ^ '^^i channel the z/^ disappearance is partially compensated 
by the conversion of (initially present in the beam) into z/^. The "suppression" of the 
muon disappearance is also influenced by matter effects, that are responsible of the left- 
right asymmetry in the plots (see [0]). Also in this case, the bounds from sub-GeV and 
multi-GeV data ("S-I-Mu" and "S-|-Mb" panels) are only slightly altered by the inclusion of 
upward-going muon data ("S-|-Mu-|-U" and "S-|-Mb-|-U" panels, respectively). 

In conclusion, the upward-going muon data tend to disfavor large mixing in the z/^j ^ z/^ 
channel and, more weakly, in the z/^ ^ channel. This indication reinforces the exclusion 
of large mixing required by the negative sub-GeV results from NUSEX and Frejus. However, 
the indications for "intermediate" mixing derived by a global fit to the atmospheric data 
are not dramatically altered by the inclusion of the upward-going muon data. In particular, 
the z/^ ^ Uj. oscillation mode is still compatible with the "minimum sacrifice scenario" of 
Refs. which is thus able to accommodate all present neutrino oscillation data, with 

the only exception of the non-fiat zenith angle distribution measured in the Kamiokande 



multi-GeV experiment PU 



10 



B. Three-flavor oscillations 

Figure ^ reports the results of our analysis of upward-going muon data in the three- 
flavor space (m^, tan^0, tan^ip), shown as planar sections at representative values of m^. 
We recall that, in each (tan^-?/', tan^ 0) panel, the right side corresponds (asymptotically) to 
pure Ufj, ^ z/e oscillations, while the lower side correspond (asymptotically) to pure ^ 
oscillations. The excluded regions in Fig. ^ connect continuously these two subcases and 
show that also genuine three-flavor cases are disfavored by the available upward-going muon 
data. These cases correspond to large neutrino mixings. In particular, the threefold maximal 
mixing scenario |^ , |50| , obtained for (sin^-?/', sin^ 0) = (1, 1/2), is rather strongly disfavored 
by the upward-going muon data. 

The exact shapes of the excluded regions should be taken with some caution, since we 
have checked that they are rather sensitive to variations in the data and in their statis- 
tical treatment; therefore, new and more precise upward-going muon data from the Su- 
perKamiokande experiment might alter significantly the bounds shown in Fig. ^. However, 
the following qualitative features seem to be relatively stable: 1) weaker bounds in the 

^ channel; and 2) exclusion of large neutrino mixing and, in particular, of three-fold 
maximal mixing. 



C. Combination of all atmospheric data 

Figure |] shows the results of a combined analysis of all the available sub-GeV, multi-GeV 
(binned), and upward-going muon data. As for Fig. |]we recall that, because of the quality of 
the available upward-going muon data, these fits might change significantly with new, more 
accurate data. Nevertheless, some interesting, qualitative features are seen to emerge: 1) in 
general, Vf. ^ oscillations tend to give a better fit than z/^ ^ 1/^ oscillations, as a result of 



the combined preference of multi-GeV data |31] and of weaker bounds from upward-going 



muon data for z/g transitions; and 2) the overall oscillation fit is reasonably good even 
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for as small as ~ 5 x 10"'* eV^. These indications suggest that oscillation signals might 
be seen better in long-baseline reactor experiments, that can probe the u^, disappearance 
mode down to ~ 10"'^ eV^. If no oscillation signal is found in such experiments, large 



horizontal stripes will be excluded in the panels of Figs. and R [32 



In conclusion, both the positive and the negative indications for flavor oscillations com- 
ing from all the available atmospheric neutrino data (including upward-going muons) can 
be made compatible, within the uncertainties, at intermediate values of the neutrino mixing 
angles. In fact, although the positive indications for u oscillations favor large mixing, the 
negative results from upward-going muon experiments (as well as from some sub-GeV exper- 
iments like Frejus and NUSEX) , tend to drag the fit towards small neutrino mixings and to 
stabilize it at intermediate values. Of course, the overall compatibility of the different data 
sets might be no longer guaranteed if the total uncertainties were smaller. Therefore, it is 
very important that the magnitude of the experimental systematic errors in the muon flux 
distribution (see Sec. II B) is clarified in the running underground experiments (MACRO, 
Baksan, and SuperKamiokande). 



IV. SUMMARY AND CONCLUSIONS 

We have performed an updated analysis of the available upward-going muon data (which 
disfavor neutrino oscillations), and a comparison with the atmospheric sub-GeV and multi- 
GeV data (which globally favor neutrino oscillations). We have worked out the oscillation 
constraints placed by the observed upward-going muon zenithal distribution, and identified 
the zones of the two- and three-fiavor mass-mixing parameter space where all the atmo- 
spheric data can be reconciled. In general, the z/^ ^ Ue oscillation channel provides a better 
fit than the z/^ ^ u^- channel, thus privileging searches with long-baseline reactors. Large 
neutrino mixing, including threefold maximal mixing, is generally disfavored. The "mini- 
mum sacrifice" fit of Refs. P^psj] , that discards, in part, the atmospheric multi-GeV data, 
is not significantly altered by the inclusion of the upward-going muon information. 
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TABLES 

TABLE L Detector characteristics, experimental data, and theoretical expectations for up- 
ward-going muons. Theoretical expectations refer to our calculations using Bartol neutrino fluxes 



|, MRS(G) structure functions [^,35], and Lohmann-Kopp-Voss muon energy losses in standard 
rock |Ho|. See the text for details and comments. 



Experiment 



Kamiokande |19| 



Baksan |22,21 



MACRO [24 



IMB-1,2 [SSl"^ SuperKamioka [ 2f 1 



Depth (m.w.e.) 
Live time (yr) 
Technique 
Dimensions (m^) 
Threshold 
No. of observed ^'s 



2700 
- 7 
Cherenkov 
7r7.82 X 16.1 
> 3 GeV 
364 



850 
11.94 
scintillator 
17 X 17 X 11 
> 1 GeV 
558 



3700 
3.06 

scint.+tracking 
12 X 77 X 9 
> 1 GeV 
255 



1570 

2.53 
Cherenkov 
18 X 17 X 22.5 
E^ > 1.8 GeV 

430 



2700 
0.63 
Cherenkov 
7rl6.92 X 36.2 
E^>6 GeV 
267 



muon flux'' 



muon flux'' 



muon flux" 



muon flux" 



muon flux" 



bin 


expt.ito-stat 


thoo. 


Oxpt.±crstat 


theo. 


expt.±(Tatat 


theo. 


expt.±CTstat 


theo. 


expt.±(Tstat 


theo. 


[-1, -0.9] 


1.44 ±0.27 


1.54 


2.01 ±0.22 


2.04 


1.10 ±0.22 


2.04 


1.72 ±0.29 


1.77 


0.86 ±0.25 


1.24 


[-0.9, -0.8] 


0.89 ± 0.20 


1.62 


2.54 ±0.27 


2.13 


1.03 ±0.23 


2.13 


1.03 ±0.23 


1.85 


1.35 ±0.31 


1.30 


[-0.8, -0.7] 


1.54 ±0.28 


1.70 


2.73 ±0.27 


2.24 


2.22 ±0.35 


2.24 


1.78 ±0.30 


1.95 


0.98 ±0.25 


1.37 


[-0.7, -0.6] 


1.44 ±0.27 


1.81 


2.43 ± 0.30 


2.37 


3.23 ± 0.46 


2.37 


1.49 ±0.28 


2.07 


1.78 ±0.34 


1.47 


[-0.6, -0.5] 


1.25 ± 0.25 


1.94 


1.82 ±0.30 


2.53 


1.54 ±0.34 


2.53 


1.97 ±0.30 


2.21 


1.26 ±0.28 


1.57 


[-0.5, -0.4] 


1.41 ± 0.27 


2.11 


2.70 ±0.34 


2.73 


2.68 ±0.51 


2.73 


2.16 ±0.33 


2.40 


1.63 ±0.31 


1.72 


[-0.4, -0.3] 


2.24 ±0.35 


2.33 


2.16 ±0.34 


3.00 


3.12 ±0.64 


3.00 


2.60 ±0.36 


2.65 


1.32 ±0.29 


1.91 


[-0.3, -0.2] 


2.69 ± 0.38 


2.59 


4.60 ±0.57 


3.31 


4.57 ±0.94 


3.31 


2.86 ±3.93 


2.93 


1.66 ±0.31 


2.13 


[-0.2, -0.1] 


3.94 ±0.46 


2.93 


5.50 ±0.80 


3.72 


4.01 ± 1.21 


3.72 


3.93 ± 0.49 


3.30 


2.71 ±0.43 


2.43 


[-0.1, 0] 


2.27 ±0.33 


3.52 


6.00 ± 1.37 


4.39 


12.7 ±5.18 


4.39 


2.71 ±0.49 


3.93 


4.06 ±0.52 


2.93 


Total muon flux" 


1.91 ±0.10'' 


2.21 


3.25 ±0.19'^ 


2.85 


3.62 ±0.55'= 


2.85 


2.22 ± 0.11*" 


2.51 


1.76 ±0.10'" 


1.81 



^Units: 10 ^"^ cm ^ s ^ sr ^ . Total flux (bottom row) = sum of partial fluxes in each cos9-bin times the bin width. 

Experimental systematic errors not reported in jl9[|5[||,§| . 
•= Experimental systematic error of the total flux estimated to be 8% 

''lMB-1,2,3 | p6) (3.6 yr live time) did not publish the muon flux distribution. The 617 observed muons were instead divided 
in "stopping" (85 events, 1.4 < Bp < 2.5 GeV) and "passing" (532 events, E^ > 2.5 GeV). The total muon flux reported in 
|26[ is 2.22 ± 0.11 X 10^^^ cm"-^ s~^ sr" ^ (stat. error only). Our theoretical estimates: 102 (stopping) + 645 (passing) muons; 
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FIGURES 

FIG. 1. Compilation of upward-going muon data and expectations for different experiments: 
Kamiokande [||, Baksan [PJH, MACRO ||], SuperKamiokande [|| IMB-1,2 [||, andIMB-1,2,3 



|26|. Dots with error bars: Experimental data with la statistical uncertainties. Solid histograms: 
our calculations. Dashed histograms: Experimental simulations. 9 is the zenith angle. See also 
Table I. 

FIG. 2. Comparison of all upward-going muon spectra, rescaled to a common threshold of 1 
GeV. (a) Separated data, (b) Combined data (as used in our analysis). See the text for details. 

FIG. 3. Results of the oscillation fit in the two-flavor oscillation channel i/^ <-> shown as 
C.L. contours in the {tan'^'tp, m?) plane. The allowed regions are marked by stars. 

FIG. 4. Results of the oscillation fit in the two-flavor oscillation channel <-> fe, shown as 
C.L. contours in the (tan^ rm?) plane. The allowed regions are marked by stars. 

FIG. 5. Bounds placed by all upward-going muon data in the three-flavor oscillation space 
(m^, tan^ (f), i&v? ip). The allowed regions are marked by stars. 

FIG. 6. As in Fig. |5|, but with the addition of the atmospheric sub-GeV (S) and multi-GeV 
binned (Mb) data (analyzed as in j31|). 



19 



Compilation of upward-going nnuon data 
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FIG. 1. Compilation of upward-going muon data and expectations for different experiments: 
Kamiokande @, Baksan [^Q, MACRO @, SuperKamiokande [||] IMB-1, 2 andIMB-1,2,3 
|26|. Dots with error bars: Experimental data with la statistical uncertainties. Solid histograms: 
our calculations. Dashed histograms: Experimental simulations. 6 is the zenith angle. See also 
Table I. 
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FIG. 2. Comparison of all upward-going muon spectra, rescaled to a common threshold of 1 
GeV. (a) Separated data, (b) Combined data (as used in our analysis). See the text for details. 
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z^^Oi/^ oscillations 




FIG. 3. Results of the oscillation fit in the two- flavor oscillation channel ^ Ut, shown as 
C.L. contours in the (tan^ ■0, m?) plane. The allowed regions are marked by stars. 
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i^^Oz/e oscillations 




FIG. 4. Results of the oscillation fit in the two-flavor oscillation channel <^ fg, shown as C.L. 
contours in the (tan^ (j), w?) plane. The allowed regions are marked by stars. 



23 



3z/ oscillations, upward— going jj. 




FIG. 5. Bounds placed by all upward-going muon data in the three-flavor oscillation space 
(m^, iw? (f), tan^V)- The allowed regions are marked by stars. 
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3iy oscillations, U + S + Mb data 
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FIG. 6. As in Fig. ^, but with the addition of the atmospheric sub-GeV (S) and multi-GeV 
binned (Mb) data (analyzed as in ||3l|). 
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